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We propose a theoretical framework based on SU{3) coherent states as a convenient tool to 
describe the collective state of a Bose-Einstein condensate of spin 1 atoms at thermal equilibrium. 
We work within the single-mode approximation, which assumes that all atoms condense in the same 
spatial mode. In this system, the magnetization ruz is conserved to a very good approximation. 
This conservation law is included by introducing a prior distribution for and constructing a 
generalized statistical ensemble that preserves its first moments. In the limit of large particle 
numbers, we construct the partition function at thermal equilibrium and use it to compute various 
quantities of experimental interest, such as the probability distribution function and moments of 
the population in each Zeeman state. When N is large but finite (as in typical experiments, where 
N ~ 10^ — 10®), we find that fluctuations of the collective spin can be important. 

PACS numbers: 67.85.Fg,67.10.Fj 


INTRODUCTION 


Coherent states (CS) are an essential tool of modern 
physics. The original (or “standard”) coherent states of 
a harmonic oscillator are quasi-classical wavepackets fol¬ 
lowing closely the classical oscillating trajectories with 
minimal uncertainty in their position and momentum. 
Mathematically, they are obtained by acting with a dis¬ 
placement operator £>(7) on the ground state, with 7 a 
complex number indexing the CS. This definition can be 
generalized to other systems, in particular if the Hilbert 
space H of the physical system under consideration is an 
irreducible representation space for a particular Lie group 
G = {(5(7)}-)- indexed by a continuous label 7. Following 
[1, 2], generalized CS are then obtained by acting with 
all elements of the group on some reference state |ref) in 
PL. The operators (5(7) generalize the displacement op¬ 
erators introduced above. The group structure ensures 
that the family of generalized CS generated in this way 
is closed. When the group is SU{2) and the representa¬ 
tion space the Fock space with N particles in two modes 
(which describes, for instance, an ensemble of N two- 
level atoms or spin 1/2 particles), this construction leads 
to the well-known SU(2) CS, sometimes simply called 
“spin coherent states”. Spin CS are used extensively in 
fields as diverse as quantum optics [3] or magnetism [4]. 
Similarly, one can introduce SU{3) CS to describe the 
collective state of N three-level atoms [5-7] or spin 1 lat¬ 
tice models of magnetic materials (e.g. [8, 9]). 

In this article, we apply the SU{3) CS formalism to 
the study of the equilibrium properties of a spin 1 Bose- 
Einstein condensate at finite temperature, where the 
three modes are the three Zeeman states with magnetic 
quantum numbers m = 0, ±1 along a given quantiza¬ 
tion axis z. We assume the validity of the single-mode 
approximation (SMA), which considers that all atoms 
condense in the same spatial mode [10-18]. Within the 


SMA, the resulting formalism allows one to describe an 
arbitrary collective spin state of N bosons in the same 
spatial mode, and is sufficiently simple to yield explicit 
analytic predictions that can be used to analyze experi¬ 
mental results. We discuss in detail the expected thermal 
equilibrium for the situation of current experimental rel¬ 
evance, where N is large and where the magnetization of 
the system is conserved [19]. We illustrate the method 
by calculating the first moments of the m = 0 popula¬ 
tion as well as its probability density function in various 
cases of experimental interest. In the thermodynamic 
limit A —>■ 00, the theory reduces to the usual mean-field 
treatment [11, 12, 20]. When N is large but finite (as in 
typical experiments, where N ~ 10^ — 10®), we find that 
fluctuations of the collective spin can be pronounced and 
experimentally observable. 


SU{3) COHERENT STATES 


In order to establish notations, we start by consid¬ 
ering a single spin 1 particle. An arbitrary pure state 
1C) = U\m = -1-1) can be obtained by applying a SU{3) 
transformation U to the maximally polarized state. We 
find convenient to express jC) as 
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/^ 1-no+mz 


l—np — mz 


( 1 ) 


The global phase is chosen such that the m = 0 com¬ 
ponent is real. Physically, uq G [0,1] corresponds to the 
reduced population in m = 0 and rriz G [—1,1] to the 
reduced longitudinal magnetization, both normalized to 
the total particle number (0 < 1 — no — \mz\ < !)■ We 
have also introduced two relative phases, 0 = £+1 + 
(/)_i — 2(/o G [0,27r] and a = (j)^i — (p-i G [0,47r], with 
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the argument of the component Cm- The phase a can 
be related to the orientation in the x — y plane of the 
average transverse spin vector sj_ = (Sxje^ + {sy)ey, and 
the phase 0 to the magnitude of s_l, 

= 2no — no + \/(l — — ml cos 0^ . (2) 

We now turn to the case of N bosons. The Fock space 
for N bosons with three possible internal states remains 
a representation space for SU{3), and we can generate a 
family of states \C^) by a SU{3) transformation acting 
on the maximally polarized state |iV : m = +1) [7]. The 
states IC^)) explicitly given by 

0 = ^ «-»’)" |e>). (3) 

in second quantized notation, describe an ensemble of N 
bosons condensing simultaneously in the same spin state 
IC). Here a = (a+i, oq, O-i)^ is a vector notation for the 
annihilation operators in each Zeeman substate, and |0) 
is the vacuum state. 

The SU{3) CS |C^) are generalized CS in the sense de¬ 
scribed in the Introduction. They provide a resolution of 
the identity operator [1, 2, 5-7, 21], f |C^)(C^I = 
1 , with the measure dC = {N + 1){N + 2)/8n^ x 
dnodrrizdQda. This implies that SU{i) CS form a ba¬ 
sis of Fock space, which is overcomplete since CS are not 
orthogonal to each other. Operators acting in Fock space 
can be represented as diagonal operators in the CS basis 
[1, 2]. 

Classical limit N ^ oo 

The main interest of using the basis of coherent states 
is the simplicity of the resulting theory in the large N 
limit, as discussed in details by Yaffe [21] (see also [15] 
for a discussion focusing on Bose-Einstein condensates). 
In this limit, the scalar product (C^IC^’^) = iC ' C)^ 
becomes very peaked around C — C due to the large N 
power. This allows one to perform the approximation 

( 4 ) 

with A a normalization coefficient. The function 6^^'^ is 
normalized to unity with the measure dC and vanishes 
very quickly (on a typical scale ^ 1/V^ in each gen¬ 
eralized coordinate) when moves away from As a 
result, tends to a Dirac distribution (5(C — C) when 
N ^ oo. The normalization of the CS implies the re¬ 
lation / dC' — CO]^ = which can be used to 

evaluate the normalization coefficient. 

In the large N limit, the quasi-orthogonality between 
two CS expressed in Eq. (4) greatly simplifies the com¬ 
putation of expectation values. In the CS states basis, a 


fc—body normally ordered operator O with k Ad N can 
be approximated as 

(C^IOIC'"^) « (C'^IOIC'^)A ^W(C - C') (5) 

As a result, the expectation value in a CS of a prod¬ 
uct of two (few-body) operators obeys a simple rule, 
(C’^IO.PIC^) ~ (C^|0|C'^)(C^I-P|C^)- This property al¬ 
lows one to compute thermodynamic averages in a ther¬ 
mal state described by a density operator p using the 
intuitively appealing formula, 

(o) ^ jdc X (c'^ipk^). (6) 

Effectively, (C'^IpIC^) plays the role of a classical distri¬ 
bution function in ^ space. For the canonical ensemble, 
for instance, we have 

(7) 

with Z = J the partition function and 

with /3 = IfksT the inverse temperature. 

In this article, we consider particle numbers N which 
are large but finite, so that the orthogonality relation 
between two coherent states holds only approximately. 
In the following, we use systematically the large N limit, 
which should thus be understood as the dominant power 
of A^ in a 1/A^ expansion. 

APPLICATION TO THE STATISTICAL 
MECHANICS OF A SPIN 1 EEC WITH 
CONSTRAINED MAGNETIZATION 

We now apply this formalism to the description of the 
low temperature properties of a trapped gas of ultracold 
spin 1 bosons. As indicated in the Introduction, we as¬ 
sume that all atoms occupy the same spatial mode (^(r) 
(but not necessarily the same spin state) [16]. The Hamil¬ 
tonian describing the spin dynamics is then [10] 

H = - qNo, ( 8 ) 

where Ug is the spin interaction energy per atom, g > 0 is 
the quadratic Zeeman energy, S is the total spin operator, 
and Nm is the number operator for the Zeeman state m. 

In most experiments with spinor gases {e.g. [22-28]), 
the system is prepared with a prior distribution that de¬ 
pends on the particular experimental sequence, noted pm 
{M being the eigenvalues of iSj). Typically, pm is peaked 
around the average value = Nfriz. Because Sz com¬ 
mutes with the interaction Hamiltonian, the prior distri¬ 
bution is essentially preserved by binary collisions driving 
the gas towards kinetic equilibrium. In order to account 
for the experimental situation, we use a generalized sta¬ 
tistical ensemble, where the energy and the populations 
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Pm are conserved on average. Maximizing the entropy 
following the standard Gibbs procedure leads to the den¬ 
sity matrix 

M=N 

^ Pm: (9) 

M^-N 

where Pm is the projector on the subspace M and where 
the pm's. are Lagrange multipliers introduced to enforce 
the conservation of the probabilities pm = {Pm)- 

Although this procedure would be the most rigorous 
one, it leads to a rather complicated formalism, where 
N constants of motion are required to describe the en¬ 
semble. Instead of constraining the full distribution, we 
choose in the following to constrain only the first two 
moments and Am\ = We expect that the 

differences from the more rigorous formalism will not be 
significant as long as only few-body observables are com¬ 
puted. Constraining the first two moments of leads 
to a density matrix 

P = 1 ^ 1 ( 10 ) 

where A 1/2 are two Lagrange multipliers. 

In the large N limit, the partition function is de¬ 
termined by the free energy K = Using 

Eqs. (2,8) and the properties of CS, we rewrite K as 

I3K = Y “ ^^0 ( 11 ) 

-I- /3'no (l-no + V^l^-no)^ — ml cos( 0 )^ . 

We have introduced two dimensionless parameters 

77 = NPq, /3' = NpUs, (12) 

as well as two new Lagrange multipliers /3' = 13' + 
2 N‘^X 2 and m* = —NXi/P'^ determined by the 
two constraints friz = ^ J d^mze~^^ and Am^ — 
/ dC (ruz — TOz)^ . The parameter m* is approxi¬ 
mately equal to the average magnetization (but not ex¬ 
actly, unless friz = 0 )- In the natural energy unit of 
NUs-, the parameter /3(, can be interpreted as the inverse 
of a longitudinal pseudo-temperature characterizing the 
fluctuations of Sz in the prior distribution defining our 
generalized ensemble. In comparison, the inverse temper¬ 
ature (3' characterizes the fluctuations of the transverse 
components Sx,Sy. A purely thermal prior distribution 
of mz is characterized by /?' = (3', a narrow prior distri¬ 
bution by /3' >> (3' and a broad one by j3'^ j3'. The two 

dimensionless parameters 77 and (3' allow us to specify 
the thermodynamic state completely given the two con¬ 
straints on friz and Amz- We emphasize that 77 and /3' 
are both proportional to the total atom number N^ which 
reflects the fact that we are dealing with fluctuations of 
a collective variable. 


In principle, the set of equations above can be used to 
characterize the collective thermodynamic state for any 
values of the parameters, friz, Amz, p, P'. In the follow¬ 
ing, we will illustrate the usefulness of the SU{3) CS 
approach in particular regimes, where analytical results 
can be obtained : the mean field regime (where differ¬ 
ences from the T = 0 mean field theory are small), the 
regime of small hq for antiferromagnetic interactions, and 
the regime of strong spin fluctuations when friz = 0. For 
simplicity, we restrict ourselves to the experimentally rel¬ 
evant case where the distribution of fulfills <g; 1 
{friz can be chosen arbitrarily). To reduce the number of 
varying parameters, we consider a system at a fixed tem¬ 
perature T, atom number N and interaction strength Us 
(so that f3' is constant). We vary the quadratic Zeeman 
energy (or equivalently 77 ) and the average magnetization. 


MEAN-FIELD REGIME AT T / 0 

The thermodynamics is controlled by the behavior of 
the free energy K. For all choices of the parameters friz 
and q, K has a well-defined minimum which depends on 
the sign of the spin-exchange interaction Us- At T = 0 
the atoms condense in that minimum [20], which deter¬ 
mines the phase diagram [20, 25, 28]. For Us < Q (fer¬ 
romagnetic interactions), the minimum is obtained for 
0 = 0 and a certain value tiq > 0 which maximizes the 
magnitude of the transverse spin. For C/g > 0 (antifer¬ 
romagnetic interactions), the minimum is obtained for 
0 = TT and the value tiq which minimizes S'_l . Differently 
from the ferromagnetic case, tt-q is zero until a critical 

value qc = C/s ^1 — \/l — above which it becomes 
positive. 

One normally expects that for sufficiently low temper¬ 
atures, the system only explores the vicinity of the min¬ 
imum. The free energy then equals its value K* at the 
minimum plus small additional terms, corresponding to 
Gaussian fluctuations around the minimum with a typi¬ 
cal spread Atiq, A0 ^ 1/13' (up to coefficients depending 
on q/Us and friz)- This describes well the case of ferro¬ 
magnetic systems, where the hnite T solution is always 
close to the zero temperature one as shown in Figure 1. 
The small differences are due to the combined effects of 
fluctuations and of the spread of mz, which are included 
in the finite T calculation but not in the T = 0 one. 
The Gaussian expansion is valid provided f3' » 1, or 
ksT NUs- For typical experimental values in [25, 28], 
N ~ 10^, Us/kB ~ 2 nK and T ~ 100 nK, /3' - 200 is in¬ 
deed large. This leads to results for the thermodynamic 
observables essentially identical to the ones obtained at 
zero temperatures, up to small corrections of magnitude 

~ 1//3'. 

Antiferromagnetic systems behave differently (see 
Fig. 2). For q < qc, the value of uq is not zero and fluctu- 




4 




FIG. 1: (a): Moments of no versus rj in the ferromagnetic case 
(Us < 0) for magnetizations =0.1 (red lines and symbols) 
and fhz = 0.3 (black lines and symbols). The solid lines show 
the average populations, the dashed lines the standard devi¬ 
ations, and the small dots the expected behavior at T = 0. 
(b): magnitnde of the transverse spin. The graphic conven¬ 
tions are the same as in (a). The values P' = NpUs = 200 
and Arriz = 0.02 were used for all plots. 

ations are comparable to the mean value. This strongly 
differs from the conclusion drawn from the T = 0 theory. 
Both effects become larger when ruz —>■ 0. 

ANTIFERROMAGNETIC SYSTEMS WITH 
SMALL no 

In order to understand the regime with q < and 
T ^ 0 better, we assume that uq remains small (which re¬ 
quires a finite magnetization rUz, see below) and expand 
the free energy around 0 = tt and no = 0. This gives (af¬ 
ter integration over and 0) the partition function as 
Z' oc /p dno with an effective free energy 

/3Aeff = Bnl-A(^^-l^no + 0 (ng) , (13) 

with A = Nj3qc and B = Ay/l — fn^. For q > qc, the 
free energy has a minimum for Uq « A{q/qc — 1)/2B > 0, 
where the first derivative of Keg vanishes: We retrieve 
the mean field regime. 

The equilibrium population in m = 0 is given by 
nd|g=o « IM for 9 = 0 and by n^\q=q, « ai(l - 




FIG. 2: (a): Moments of no versus rj in the antiferromag¬ 
netic case {Us > 0) for magnetizations rhz = 0.1 (red lines 
and symbols) and rhz = 0.3 (black lines and symbols). The 
solid lines show the average populations, the dashed lines the 
standard deviations, and the small dots the expected behav¬ 
ior at T = 0. The larger white dot shows the analytical limit 
for q = 5c, no|q=9c « 0.34(1 — for rhz = 0.3. (b): 

magnitude of the transverse spin. The graphic conventions 
are the same as in (a). The values P' = NpUs = 200 and 
AtUz = 0.02 were used for all plots. 


!A^/"^ (oi « 0.34) for q = q^. For high rUz, we 
have A « /3' ^ 1, and correspondingly small population 
in m = 0. For small mz, we have A « /3'm^/2 : the 
population in m = 0 is small only if tUz ^ Ij'/P'. Note 
that when this is not fulfilled {A ~ 1), the expansion 
in Eq. (13) is not valid. The average transverse spin per 
atom also becomes finite at T 0. In the limit where ng 
remains small, we find from Eq. (2) 


si - + O (nl ngO^) . (14) 


In the limit g —?► 0, we find s\ « 2//3' independent of 
m^. For the typical values of P' = 200 given above, s5_ 
reaches a few percent below q^ (see Figure 2b), which is 
experimentally measurable. 
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SPIN FRAGMENTATION FOR 
ANTIFERROMAGNETIC INTERACTIONS AND 
SMALL MAGNETIZATIONS 

For antiferromagnetic interactions {Us > 0), a spe¬ 
cial situation occurs near niz = 0 where the critical 
vanishes. In this regime, arbitrary large fluctuations of 
Uq are possible, which makes the approximation used in 
the previous paragraph invalid. Going back to the spin- 
dependent free energy Eq. (11), we expand around 0 — tt 
and rriz = 0, 



no 


1 ?n^ (0 — 7rl^ 

/3A«-^+/3'no(l-no)^^^^-r?no, (15) 

2 ( 7 ^ 2 


with 


cr = 


1 - no 


P'no + I3{{1- no)' 


(16) 


The overall minimum, determined by the quadratic Zee- 
man energy [last term in Eq. (15)], is at Uq = 1, 0 = tt 
and niz = 0. The curvature near the minimum along 
the no direction vanishes. In this case, the condition 
/I' ^ 1 is not sufficient to ensure the distribution is 
peaked around the mean-field solution: One must also 
have f] ^ 1. When <7 = 0, this is never fulfilled 
and instead of a single minimum, one finds instead a 
family of degenerate minima, corresponding to the so- 
called polar (or “spin-nematic”) states with vanishing 
spin [11, 12, 15, 29]. Because of the broad distribution, 
the system displays in the limit rj 1 large fluctuations 
in the individual populations no,n+i,n_i. As discussed 
in Refs. [13-15, 17, 29, 30], this is a signature for frag¬ 
mentation of the condensate which can occupy any of 
the quasi-degenerate states or an arbitrary superposition 
of them. We stress again that these fluctuations are a 
mesoscopic effect, and disappear in the thermodynamic 
limit where they are confined to a vanishingly small win¬ 
dow around ij = 0. Eor /?' = /3', we recover the previous 
results obtained at finite temperatures [29]. The present 
theory is able to go further by accounting for the most 
general situation where /3(, /?'. 

After integration over 0,mz, we obtain the (unnor¬ 
malized) marginal distribution function of uq in a simple 
form, 


pTino 

P{no) oc (17) 

^/noiP'no + ^'(1 - no)) 

Other marginal distributions {e.g. for m^) could be ob¬ 
tained in a similar way. Eor q = 0 (see Eigure 3), the 
distribution of uq changes from an asymmetric character¬ 
istic square-root shape, P{no) c>c \/y/no [31], to a sym¬ 
metric shape P{no) oc l/-y/no(l — Uq) when /?' changes 
from 1 (unconstrained prior distribution of niz) to -l-oo 
(narrow prior distribution of m^). Eor small ry <|; 1, the 


FIG. 3: Probability density P{no) for the reduced population 
in m = 0 (the densities are not normalized). The red solid 
line shows the distribution for = P' and the blue dashed 
line for /3( = 100/3k 



FIG. 4: Approximate theory in the case rfiz = 0 and q = 
Average (solid line) and standard deviation (dashed line) of 
the population no. The vertical dotted line marks the uncon¬ 
strained case [P'z = P'). 


distribution P{no) is always broad, so that the qualita¬ 
tive conclusions about spin fragmentation and large pop¬ 
ulation fluctuations are unchanged. 

To discuss the influence of the prior distribution of , 
we set <7 = 0 and plot in Figure 4 the first two moments 
of no versus /3'. For /?' = 1, which corresponds to the 
situation without constraint where ug = 1/3, we find 
Ang « 0.30 and Aniz = 2/{PP'). With increasing /3', 
the prior distribution of ruz becomes narrower: ng goes 
from 1/3 to 1/2 and Ang increases slightly [29]. Con¬ 
versely, /3' < 1 corresponds to a prior distribution of niz 
broader than the one without constraint : the average ng 
decreases below 1/3 (the fluctuations of ng decrease as 
well). 


CONGLUSION 


In conclusion, we have proposed a description of the 
collective equilibrium state of a spin 1 Bose-Einstein 
condensate based on SU{3) coherent states. Using this 
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formalism, the conservation of the magnetization is ac¬ 
counted for by introducing a prior distribution for 
and constructing a generalized statistical ensemble that 
preserves its first moments. We have computed moments 
of various quantities, and their probability distribution 
function (for example, for the population in the m = 0 
Zeeman substate), that can be directly compared to ex¬ 
periments. Going beyond thermodynamics as studied in 
this paper, we expect, in analogy with what has been 
done with SU{2) coherent states, that the SU{3) co¬ 
herent state formalism can be used to study the collec¬ 
tive dynamics leading to spin oscillations [32-35] or spin- 
nematic squeezing [36-39]. Another interesting direction 
is to extend the formalism to larger groups SU(N) with 
> 3 to describe condensates with higher spin. This 
is relevant for instance for experiments with Chromium 
atoms with spin 3 [40, 41], where spin exchange and mag¬ 
netic dipole-dipole interactions both play an important 
role. 

We thank the members of LKB for stimulating dis¬ 
cussions. We acknowledge support from IFRAF, from 
DARPA (OLE program), from the Hamburg Center for 
Ultrafast Imaging and from the ERC (Synergy grant 
UQUAM). 
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